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Abstract 

Magnetic nanoparticles are useful biological probes as well as ther¬ 
apeutic agents. There have been several approaches used to model 
nanoparticle magnetization dynamics for both Brownian as well as Neel 
rotation. The magnetizations are often of interest and can be com¬ 
pared with experimental results. Here we summarize these approaches 
including the Stoner-Wohlfarth approach, and stochastic approaches 
including thermal fluctuations. Non-equilibrium related temperature 
effects can be described by a distribution function approach (Fokker- 
Planck equation) or a stochastic differential equation (Langevin equa¬ 
tion). Approximate models in several regimes can be derived from 
these general approaches to simplify implementation. 


1 Introduction and nanoparticle applications 

Magnetic nanoparticles (MNPs) are useful in many biophysical and medical 
applications because they can be remotely controlled and monitored with 
magnetic fields. For these applications, it is important to model the dynam¬ 
ics of the particles; and many approaches have been used previously to do so. 
In this case ‘modeling’ means simulating the magnetizations of the particles 
over time to study how different variables (e.g., field strength, anisotropy) 
affect their dynamics. In this review, we outline the major advances in MNP 
modeling including state-of-the-art techniques requiring numerical simula¬ 
tions and accounting for many experimentally verified phenomena. We also 
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include approximate theories that admit analytical solutions, formulations 
that are easier to implement and other expressions that are more transpar¬ 
ent. 


1.1 Sensing and imaging 

Magnetic nanoparticles can be used as sensors to detect properties of their 
local micro-environments when rotations couple MNP dynamics to the en¬ 
vironmental parameters [H 121 |3j. For example, the concentration of specific 
molecules can be detected because their dynamics can be noticeably differ¬ 
ent when bound and unboundjl]. This is depicted in Fig.[^a). Nanoparticle 
sensing is sensitive (down to lOOpM of analyte and nanogram amounts of 
iron) and could potentially be used in vivo to detect changing concentra¬ 
tions over time. Aside from molecular sensing, MNPs have been used as 
probes to measure local temperatures [5], viscosities [6], and local environ¬ 
ment rigidity [7]. MNPs have been used as MRI contrast agentsjH], and the 
developing technology of magnetic particle imaging (MPI). MPI uses the 
particles themselves as high-contrast imaging agents which hold the possi¬ 
bility for time-domain tracers [9l [TOl [TT]. 

1.2 MNP hyperthermia 

Magnetic nanoparticle (MNP) hyperthermia is considered a potentially use¬ 
ful addition to current cancer treatment modalities [121 US]. Nanoparticles 
can be directly injected[Tl| or targeted biologically [15] to achieve specific lo¬ 
calized therapy. Furthermore, it appears cells ingest the MNPs, so cellular 
distribution may be possible [16]. If oscillating magnetic fields are applied 
to the sample, diverse physical mechanisms prevent the particles from fol¬ 
lowing the field exactly, and thus energy is dissipated locally. This energy 
can be thought of as heat, and it has been shown that cytotoxic heating is 
possible [T7[ ITS]. The principle of MNP hyperthermia can be conceptualized 
as in Fig. Bb). It is important to design efficient heating agents and much 
work has been done to characterize MNPs and improve properties such as 
magnetic moment magnitude and anisotropy energv fl^ [T7]. 
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Figure 1: Two applications for magnetic nanoparticles (a) sensing specific 
‘analyte’ biomolecules, and (b) locally heating through dissipative losses 
from oscillating particles. 

2 Nanoparticle rotations and relaxation times 

There are two mechanisms for nanoparticles to rotate their magnetic mo¬ 
ment /i. This is visualized in Fig. The whole particle can rotate, as in 
so-called Brownian rotation, named for Robert Brown. The moment can ro¬ 
tate internally due to restructuring of electronic states. This so-called Neel 
rotation is named for Louis Neel, who first described the phenomenon while 
studying magnetic remanence in geological samples |20|. The timescale for 
a perturbed system to return to equilibrium is called the ‘relaxation time’. 
This expression is different for each mechanism and depends on the param¬ 
eters of interest. 


2.1 Brownian relaxation time 


The Brownian relaxation time, sometimes called the Einstein relaxation 
time m based on Einstein’s original work to derive the timescale from Brow¬ 
nian motion arguments is written 


tb 


3r]Vh 

kT 


( 1 ) 
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Figure 2: Neel and Brownian particles relax differently. The moment can 
rotate internally to align with the anisotropy axis n, or the entire particle 
can rotate. 

in terms of the suspension viscosity 77 , the hydrodynamic volume of the 
particle Vh and the thermal energy introduced with Boltzmann’s constant k 
and local temperature T. 


2.2 Neel relaxation time 


The Neel relaxation time derives from the thermal movements between two 
potential energy minima that arise from postulating a single anisotropy axis 
n (see Fig. |^. The anisotropy energy is written in terms of the anisotropy 
constant K and the magnetic core volume Vc] the unitless ratio of this 
energy to the thermal energy is defined a. Thus, the relaxation time can be 
expressed 


tn = 


KVcIkT 


Toe 


( 2 ) 


where the ‘event time’ tq is often quoted as 10 ^^s but is defined 


/i (l + a^) 

2'ykT a 


(3) 


in terms of the electron gyromagnetic ratio 7 = 1.76- lO^^Hz/T, a dimension¬ 
less magnetic damping parameter a, and the magnetic moment—calculated 
using Mg the material dependent saturation magnetization as /i = MgVc. 
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More precise forms of the Neel relaxation time including additional factors 
of a can be found in Ref. m- 


2.3 WF relaxation time 

W F Brown developed a more general characteristic time including applied 
fields [23]. Letting a 1 he approximated this time to be 

_ _^N_ 

'AVF 1 — ^2 ^ e sinh ^ 

where the applied magnetic field H is accounted for with the unitless ^ = 
fiH/kT. The ratio of anisotropic to magnetic energy is defined e — 
and the expression reduces to the equilibrium Neel time in the absence of 
an applied field. 


2.4 Relaxation time considerations 


An assumption found frequently in the literature is that the more prevalent 
relaxation mechanism is the one with the shorter relaxation time. The 
process is then approximated as parallel so the effective relaxation time 
is 


1 _ 1 1 

Feff Tw Tb 


(5) 


To the best of our knowledge there has been no theoretical work that deals 
with Neel and Brownian relaxation simultaneously, so currently this expres¬ 
sion is not derived from first principle and does not incorporate the fact 
that the processes could be coupled [SUES]. Another consideration is that 
these relaxation times only hold for equilibrium conditions, and should not 
be applied without awareness of this constraint. 


3 Anisotropic hysteresis loops 

The energy of a magnetic particle with magnetization direction m and uni¬ 
axial anisotropy direction n immersed in a magnetic field H can be writ¬ 
ten 

E = -KVc (n • m)^ - /im • H (6) 
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so that if the field is applied in the z direction, we define the angle from the 
polar to be 9 and the angle between the moment and the easy axis to be 0. 


Thus 


E{e,(p) = KVcSin^ie - 0) - fiHcose. 


(7) 


To find the minimum energy, the equilibrium state of the particle, with 
respect to the moment direction 9 we write 



( 8 ) 



zation M — cos 9 and plotting this with respect to the normalized field e we 
find that the so-called Stoner-Wohlfarth hysteresis loop emerges [26]. 

This model is useful because the area of the loop indicates the energy dissi¬ 
pated when the particle is forced to change its orientation from one poten¬ 
tial minima to another. This can then be used as a model for hyperthermic 
heating given particles with specific sizes, saturation magnetizations, and 
anisotropy constants. The weakness of the model is that it assumes no re- 
laxational timescale for the process, thus neglecting that any phenomena 
depend on the frequency of the applied field. Still, models based on the 
theory have been extended to model heating more realistically [27]. 

4 Differential equations for the magnetizations 

There exist more general methods to model the behaviors of the magnetic 
moments in various applied fields. This can be accomplished by describing 
the time dynamics of the magnetization using differential equations. Then 
the varying rotation methods as well as the specific conditions the particles 
experience in various applications can be simulated. 
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4.1 Brownian dynamics 


A Brownian particle can be modeled phenomenologically with a balance of 
torques [28]: 

dm 

T = /im X H — 67714111 X . ( 10 ) 


dt 


magnetic 


( 11 ) 


We have neglected the acceleration term because the nanoparticles have a 
very small Reynold’s number. Balancing this equation admits a differential 
equation for the magnetization which can be rewritten in terms of a unitless 
field ^ = /iH/fcT and the Brownian relaxation time tb as 

dm (m X X m 
dt 2 tb 

This magnetization equation can be enhanced by including thermal ffuctu- 
ations of the magnetization of the particles. There is no one general way to 
do this, and we will see the prescription for Neel particles is slightly differ¬ 
ent. Here we supplement the torques with the stochastic term Ts — 

The magnitude of this term depends on the relaxation time so that larger 
torques occur at higher temperatures and lower viscosities. We have also 
introduced the white noise process Xt which we define to have 

(y) = 0, (12) 

where the process is uncorrelated in each cartesian direction G x^y^z 
and is called white because its variance is a Dirac delta function in the time 
domain and thus is a constant in the frequency domain m- In a physical 
system, noise is not truly white but is a fair assumption when the stochastic 
torques have a much shorter timescale than the rotations of the particles. 
We write the stochastic differential equation, called a Langevin equation [30] 
as 

dm (m X X m At X m 
dt 2 tb y/TB 

This equation is used to model Brownian particles and is amenable to chang¬ 
ing local variables like ffuid viscosity and temperature as well as externally 
applied fields. 


4.2 Neel dynamics 


The equation for the change in the magnetization of a Neel particle is the 
phenomenological Landau-Lifshitz-Gilbert (LLG) equation. This can be de- 
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rived from the Larmor precession of a spin in a magnetic field with an added 
velocity-dependent damping term m- The normalized internal moment 
then rotates with 

dm y 

-= -^ [H X m + am x (H x m)l . (14) 

dt 1 + a^ ^ ^ ^ ^ 

The effective field H is used to include additional dynamics beyond the 
applied field [32]. For example we include an externally applied field with 
amplitude Hq and frequency cj, a field from the anisotropy axis n[33|, and 
a stochastic field h, so that 

H HoZCO^ujt-\- • n)n+ h(t) . (15) 

'---" jJi 

applied ' v-' stochastic 

anisotropy 

In this case, the fluctuations are added to the model with an additional 
stochastic field, a white noise field with zero-mean and standard deviation 
parameterized by the nanoparticle variables with h(t) = Some 

work has been attempted to describe the effects of dipole fields, which can 
also be added to the effective field [3^ or treated as a mean field|3l|. But 
because it seems the effects are actually of detriment to the effectiveness of 
hyperthermia, we do not discuss them here. 


4.3 Numerical integration of Langevin equations 

With the addition of the stochastic terms, the differential equations cannot 
be integrated as a normal Riemann integral [35]. In general we see that the 
dimensions of the stochastic terms are actually the square-root of time, so in 
fact it is more precise to instead write these equations as integral equations 
of the form 

dm = J adt -h J bdWt (16) 

where the Wt represents the Wiener process, the integral of white noise 
and also called Brownian motion—the continuous analogue to the random 
walk. Additional information on the various types of stochastic differential 
equations can be found in Ref. [33 ES]- The usual method is to integrate 
the differential equation to first or second order, using the Euler-Marayuma 
or the Heun integration scheme, respectively [37] l28]. Then, magnetization 








moments are found from successive solving 


= (17) 

i 

SO that the average magnetization is when j — 1. 

5 Distribution function approach leading to ap¬ 
proximate models 

We saw in the Langevin equation formulation that it is possible to include 
thermal effects to models of nanoparticles. This is particularly useful when 
the particles are expected to be at room temperature (as in biological and 
medical applications). The finite temperature will cause magnetization ffuc- 
tuations and thus the nanoparticles are accurately described as a distribution 
of states. Instead of repeatedly solving a stochastic differential equation, we 
can use a distribution function approach where the distribution of states 
evolves over time. 


5.1 The Fokker-Planck equation 


We can represent a single nanoparticle’s magnetization as a point on the unit 
sphere, so that the surface density of many magnetizations is determined by 
a function /(0,(/),t). As the particles rotate in space, the magnetization 
surface density changes, leading to a surface current J. The total number 
of nanoparticles is conserved, i.e., J /df] = 1, defining the surface density 
function as a probability distribution. The normalization also implies the 
continuity equation 


dt 


-V-J. 


(18) 


The probability current depends on the probability itself; this is the key to 
understanding the dynamics. At temperature T = 0 the current density J 
depends only on the velocity (v = of the points on the sphere. However, 
when the temperature is not zero, the system evolves towards equilibrium. 
A postulated phenomenological ‘diffusion’ of the distribution function [23] 
accounts for the approach to equilibrium. This is represented by a new 
term in the continuity equation proportional to a diffusion constant D and 
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the gradient of the distribution function. We write the Fokker-Planck (FP) 
equation for the distribution function as 


dt 



DV 


(19) 


From this distribution, magnetization statistics can be determined using the 
definition of the probability moments 




( 20 ) 


What remains is to define the change in the magnetization and the diffusion 
constant. The magnetization dynamics are controlled by the zero temper¬ 
ature differential equations Eq. HH and Eq. The diffusion constant is 
determined by the parameters at equilibrium conditions, when there is no 
applied field and ^ = 0. In general, solutions to the Fokker-Planck equation 
are not possible, so approximation methods are used in practice. 


5.2 Neel rotation and hyperthermia 

If the expression for the magnetization dynamics of Neel rotation is inserted, 
this FP equation can be used to describe the rotations of small particles 
that are fixed spatially. These are the expected conditions for nanoparti¬ 
cles during hyperthermia [32]. It is also possible to develop an analytical 
approximation based on the assumption that the distribution function is 
linear in the magnetization. This approximation has been used to model 
hyperthermia because the imaginary response indicates heat deposit ion [38]. 
Though these expressions are commonly used, we do not describe them here 
because recent work has shown that they are inadequate to fully model heat¬ 
ing [25l ETj . The methods are accurate but care is necessary to assure that 
they are applied in the correct range of validity: when magnetizations are 
indeed linear as when applied fields are weak, magnetic moments are small 
and/or frequencies are high [39]. 
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5.3 Brownian Fokker-Planck with cylindrically symmetric 
applied field 


Replacing the velocity of the magnetization in Eq. [T^with Eq. [T^ simplified 
slightly we have the FP for Brownian rotation 


dt 


-V- 


2 r 


-DV 


/• 


( 21 ) 


A general solution is not currently analytically possible, but if the applied 
magnetic field is in the z-direction only (i.e. ^ the FP equation can 

be simplified. 


This involves writing out all the components and unit vectors of m in spher¬ 
ical coordinates, 


m = sin 9 cos (px + sin 9 sin (^y + cos 9z 


and each unit vector is expressed e.g.. 


_ . ^ d 7 1 d . . 

X = \/x = I 9— + I sm6>cos( 

o9 sm 9 0 (p 


eventually we find 


z — m cos 9 = — sin 99 


where the second term completely cancels. Thus, we arrive at the intuitive 
result that the distribution function does not depend on the azimuthal angle 
/ = f{9^t) only. We now have the FP equation 


dt 


-V- 


1 

2 r 


sinOf -D 


df 

oe 


( 22 ) 


where by using the definition of the spherical gradient over only the polar 
angle we write 


dt 


1 d 
sin 9 d9 


sine ( -^sinOf -D^ 


(23) 


and lastly with the change of variables x = cos 9 we have — W ~ §x 
with T) = 1/2t we can write the 1-D FP equation. 



A 

dx 


(1 - 



(24) 


This equation has solutions in a Legendre polynomial expansion [40|. l34] . 
But it is of interest to demonstrate a further approximation method, the 
so-called ‘effective field’ method or ‘macroscopic relaxation equation’. 
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5.4 Effective field model for low-frequency Brownian relax¬ 
ation 

If we set the applied field to be a constant and let ^ = 0, we find the 
normalized distribution function 


fix) 


47 r sinh ^ 


(25) 


This is the same result as we get from Boltzmann statistics (remembering 
X = COS0 in terms of the polar angle because the field is applied in the z 
direction). We write the average magnetization in the direction of the field 
as M. Using this distribution function and integrating to get the first mo¬ 
ment, we find, as we should, that the Langevin function C{x) describes the 
equilibrium dynamics of an ensemble of particles in an applied field: 


M = £(^o) = coth^o - 1/^0- 


(26) 


The valuable approximation is to assume the particles are always at equilib¬ 
rium so that the distribution function retains the same functional form, but 
with a time varying field ^ ^ This ‘effective field’ approach turns out to 
be very useful to model low frequency oscillating fields where the relaxation 
time of the particles is shorter than the period of the applied field [341 [28]. 
Multiplying by x and integrating Eq. with the effective distribution and 
definitions of the probability moments leads to an equation for the second 
moment in terms of the first. We have then a differential equation for the av¬ 
erage magnetizations at “near-equilibrium” (which is notoriously a slippery 
condition): 

dM^_M / 

dt \ Ceff / 

where ^eff is found from inverting the Langevin function at every step [34]. 
An inverse Langevin function can be used as in Ref. [41] . The effective field 
equation is useful for quicker analyses because it is much less computation¬ 
ally challenging than the full stochastic models. 



6 Summary 

We have seen that there are two mechanism for a magnetic particle to 
reorient its moment. The energy changes for particles to relax internally 
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over energy barriers are also discussed in the framework of localized heat¬ 
ing as is used in magnetic nanoparticle hyperthermia. The dynamics of 
these magnetic moments undergoing each relaxation mechanism separately 
are explored. The two most general approaches (the Langevin equation, or 
stochastic differential equation approach, and the Fokker-Planck, or distri¬ 
bution function approach) are illustrated. These include full variability of 
nanoparticle parameters as well as thermal ffuctuations induced in a realistic 
setting. The methods for numerical integration are highlighted, and read¬ 
ers are directed to valuable sources. However, these equations are difficult 
and time consuming to solve, so we also introduce a macroscopic relaxation 
equation that is useful for low-frequency Brownian simulations. 
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